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ABSTRACT 

A diagnostic, non-linear balanced model is applied to a 
case study in the tropical Western Pacific on 1 and 3 March 
1965. The region includes both sides of the equator and 
eontains “the weileZe 

A discussion of the ITCZ is given and the model is dis- 
cussed along with the processing of the input data. The study 
contains the horizontal wind velocity and the thermal dis- 
tribution as input values. The non-divergent stream function 
(Y) is obtained by relaxat#n- of yyhe Porsson Hauation, vy = 
IR-VxvV2S , where the vorticity is computed from the observed 
Wind field. A comparison is made between the non-divergent 
part of the wind and the actual wind. Computer analyzed 
stream functions at several levels are shown and discussed. 

Dry adiabatic vertical velocities are obtained and com- 
pared with the ITCZ. Typical magnitudes are on the order of 
0.5 cm/sec. The forcing functions of the w equation are 
Shown and the contributions, in partitioned form, by various 


terms are illustrated and discussed at several levels. 
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lL. .ANTRONGE TIS) 

The tropical atmosphere has been the object seis nai 
tensive research effort in recent years. No single, dynam- 
ically consistant model has been developed which yields conm- 
pletely acceptable results, and although great progress has 
been made through computer application to case studies, some 
of the questions raised are still unanswered. 

Lack of conventional data severly hinders analysis in 
any tropical study and it seems unlikely that a significant 
improvement in the number of reporting stations is forth- 
coming in the near future. The alternative for those in- 
terested in the problem is to rely more heavily on "uncon- 
ventional" or satellite data, including photographs and 
radiation data, and possibly satellite tracking of constant 
level balloons. The Intertropical Convergence Zone(ITCZ) 
the Equatorial Trough, and the Intertropical Front have 
been used by various meteorologists to describe the same 
meteorological phenomena. A finer distinction could per- 


maps be Made, ebutwior thempurpose of this #aper [7TCZ is 
preferred and will be defined as a relatively narrow, 


guasi-continuous band of convective clouds encircling the 
earth with associated synoptic scale vertical velocities 
caused by the convergence of air masses originating in the 
ieee Wid Saavems Of Botm nemispheres. Uhesouner defini 
tions will be used interchangeably when referring to the 
literature. The weather associated with this zone is often 
pevere,. CONVeECLI Vermin Nature. and in the tCrepiucemexcecdcd 


only by hurricanes and typhoons. 
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A multi-level diagnostic non-linear balanced model 
is a useful tool in studying the three-dimensional struc- 
ture of the atmosphere and such a model has been success- 
fully used by Krishnamurti (1966) in tropical and mid- 
latitude studies. Charney (1963) argued that vertical 
motions in the tropics are small and may be neglected. 
Krishnamurti believes that even though small, they may 
play an important role in the dynamics of the tropics. 

He found typical values of vertical velocity on the 


1 


order of 0.1 cm sec ~ in a typical Carribbean easterly 


: in the considerably more active 


wave and 1-2 cm sec 
Indian monsoon circulation which becomes part of the 
ITCZ during the northern summer. Together with the 
horizontal winds, the vertical velocity completely 
describes the three-dimensional motion of the atmos- 
phere. 

The authors applied the model used by Krishnamurti 
and Baumhefner, (1966) with some modifications, to a 
case study in the Western Pacific in order to obtain 
vertical velocities in the ITCZ. Other considerations 
motivating the work were, to investigate the non-diver- 
gent stream function derived from the relative vorticity 
field on and near the Equator, and to define initial 
state conditions for use in a more sophisticated model, 
such as a primitive equation prediction scheme as sug- 
gested by Charney, (1962), 

A barotropic model is currently being used to fore- 


cast the nonedivergent wind field and at present is pro- 


NZ 


ducing acceptable results. Once the data problem in the 
tropics is solved, however, the baroclini ceteris: 

have to be included and a Significant improvement should be 
realized. 

History. Historically the IICZ has been a "frustrating 
and confusing phenomenon. Lacking a complete understanding 
of the major forces involved in the tropics, meteorologists 
of various schools attempted to apply their theories to the 
tropical problem. Palmer (1951) and Riehl (1954) divide the 
approaches into three broad categories: Climatological, 
Air Mass, and Perturbation. The Climatologists based their 
interpretation of the ITCZ on long term mean charts, statis- 
tical in nature, and based on a simple general circulation 
theory. The data network in the tropics may be sufficient 
for their purposes, however, they are unable to explain the 
large daily fluctuations and severe storms frequently found 
in the equatorial regions. The mean position of the equa- 
torial trough, based on isobaric analysis, is shown in 
Fig. 1. The Scandinavians, who are primarily responsible 
fer the air mass approach, applied their frontal theory to 
the tropics and the term Intertropical Front was coined. 
Indeed, the weather was of a frontal nature caused by the 
air masses of the two hemispheres converging in a low 
pressure trough, although a significant density discon- 
tinuity was not observed. The ITCZ was subsequently de- 


scribed by various authors as: 


S 


"ACW ATAee uy eysy JSeouqNnog LaeAo pue gg ATenIqeg 
pue 2G 4Snzny jo syuquom ZDT 9uq ut YsNoI4 [eyzloyenby” euy jo SUOTITSOd uBsoy °T °Stq 
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1) a single and continous front 

2) a double front 

3) part single and part double fronts 

4) discontinuous 

5) sloping toward the summer hemisphere 

6) sloping arbitrarily 

7) not silcpmunerat all 

8) it moves discontinuously 

9) it disappears at one location while 

another forms at a new location. 

It is not surprising that this theory also failed the test 
in that it could not be used as a short term forecast aid. 
No dynamic explanation could be offered for the varied 
structure and movement of the front and the theory was 
abandoned. 

Riehl (1954) describes Palmer's attempt to consol- 
idate the previous work into a unified picture in the 
Eee vurbation approachegieem,. In ar Pacific studyehe 
found that the East Central Pacific wind constancy is 
high and the equatorial trough has little variation. 
small disturbances that perturb the flow are either 
brought into the area or develop there. They then move 
downstream as waves in the Equatorial Easterlies. Con- 
vergence in the basic flow tends to amplify cyclonic 
crests and dampen anticyclonic ones. If the axis of 
symmetry for the waves is on the equator, both northern 


and southern crests of the waves represent cyclonic 


i 


circulation and will intensify due to the convergent 
flow. This may lead to what has been termed a split or 
double equatorial trough. With continued amplification, 
closed circulations appear and a west wind finally appears 
at the equator. The result of a large number of these 
circulations near the equator becomes evident in a mean 
wind analysis showing a west component in very low lati- 
tudes. These west winds are shown by Wiederanders (1961) 
to appear only below 500 mbs. Originating in the trades, 
they cross the equator, recurve and form a westerly 
asymptote of convergence. This phenomena seems apparent 
In Ghis sud, on i Maren. 

The Perturbation approach, although a signif icant 
improvement over the Climatological and Frontal approach, 
was not wholly satisfactory and many of the dynamic prob- 


lems are still with us. 
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Tl. DESCRIPTION Of arb SMe mE 

Dynamical multi-level models enable a study of ver- 
tical motions of a particular system or area. The concept 
of using balanced models for providing an initial state 
for primitive equation models was first proposed by 
Charney (1955). The balanced model initial state has the 
advantage of suppressing the unwanted gravitational in- 
stabilities that observed initial data might generate. 

In recent studies of tropical storms, balanced models 
have been proposed by Ooyama (1963), Charney (1963), 
Rosenthal (1964) and Kuo (1965). 

A multi-level diagnostic balance model has been 
demonstrated to be a powerful tool for studying the 
three-dimensional details of the non-geostrophic be- 
havior of the atmosphere during various stages of a 
storm development (Krishnamurti, Nogues and Baumhefner, 
1966). Baumhefner (1966) used a similar model, but with 
additional terms, for a case study in the Caribbean. 

The model presented here is essentially the same as that 
used by Baumhefner, but with minor modifications in cal- 
culating the non-divergent stream function and using 
actual observed temperatures, as suggested by him. 

It DMst be Kepiaienmmnd= that ethere is actually a 
two-fold purpose involved with the use of this model. 

One is diagnostic, and the second, and more far reaching, 
is that these results are to be used in the future ina 
multi-level primitive equation prediction scheme as input 


data. 
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The balance model provides a synoptic scale vertical 
motion field which is determined by numerical calculation on 
a horizontal mesh whose size is a few degrees of latitude. 
The numerical model selected to describe the three-dimensional 
structure of the tropics is commonly called the non-linear, 
or multi-level balanced model (Charney, 1962). Baroclinic 
vertical motions of a complete balanced model can be parti- 
tioned to determine the unique contributions of several for- 
cing functions. These may be looked on as several synoptic 
mechanisms producing ascending or descending motions. 

A formal derivation of this model has been done by 
Baumhefner (1966), and further explanations are contained in 
Krishnamurti (1966b). Therefore, this work will not be re- 
peated here, other than to give the final form of the omega 
equation of a general balance model and to point out the al- 
ternate method that was used to calculate the non-divergent 


stream function. 
The ®&-equation of a general balanced model is expressed 


by the following three equations for W,% , and = - 


Youws+fze= fEIWS)+ 09's (4,8) 
-23t at 3) T 55 (SV°X) 
we > 935/50  - -=— Vv" H. 
Rv*H, + f2,{ox vy} 
io 72 ae) - f 3% {vx V3.3 
~nv*{vx-'7e}-B4s & ot 4) 


ow 


7*X= oP (2) 


3 alee an i 
Voge = 7 TY S)- 95, - BF] 
+ VX-°V§, +S. V'*X 

2 
- Vo -Vig- wo 2 ¥ Wy. (3) 


These equations are solved by numerical techniques. 
The general balance W -equation, presented here, has 
liz2ehore Ings 1unct Loms: 
The forcing functions are: 
ait. Toeeo CUSt) Differential vorticity advection by the 
non-divergent part of the wind. 
Zoe a vy a(¥, Q) Laplacian of thermal advection by the 
non-divergent part of the wind. 


3. -23¢ % IS, 3H) Dif temcmutal derormavien erfect. 


2 
3, (SVX) Differential divergence effects of a 
balance model. 


f 3.9 5p LH Lot} 
De 4 ox Sy Effects of frictional stresses. 


2 
cE ee Vo Ay Effects of latent heat. 
2 
6 -& V Hs Effects of sensible heat transfer from 
P 


Water surfaces to the atmosphere. 
2 s45% Ta 
oe Tap op ¥¢ Differential vertical advection of 
Violence | Gay ¢ 


my 
9. fi jve vt? Differential turning of vortex tubes. 


s 
10. = {VX-V Sa} Differential advection of vorticity by 
the divergent part of the wind. 
2 
‘bee i 171: V9} Laplacian of thermal advection by the 


divergent part of the wind. 
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1 
alee —-Bses 


qje¢ 
ct 


Contribution by the beta term of the 
VOruLeitcy equation. 

It should be noted here that, in this form, none of the 
terms "blow up" at the equator where the Coriolis parameter 
(f) goes to zero. Instead, the forcing function just drops 
out. Zero vaiues are not found in the output, however, due 
vO Une Lormwomwmyic rolaxealLion process, 

Calculations of the input values to solve our equations 
are, for the most part, pretty svraightforward. Five point 
schemes are used for the Laplacian, the Liebmann forward ex- 
trapolation technique of relaxation is used, Jacobians are of 
the standard form and all derivatives are evaluated over dis- 


tances 24x and zay. The following equations are pertinent. 


Potential temperature: = 1 ea ep (4) 
T 3 

Dips abi.or scabs lity: =~ = 3p (5) 
oH 

Absolute vorticity: Se = a ay +f (6) 


Since an analysis of the geopotential height field is 
not easily possible in the tropics, the so-called balance 
equation (7), which is a simplification of the divergence 


equation, was not used. 
V: fVy= V* - 23 vee (7s 


If the balance equation is used, the procedure is to 
find the non-divergent stream function ( y ) from the geo- 
potential height. The stream function is then used to find 


the vorticity. 


fa) 


The procedure was reversed, somewhat, inwehis paper. 
The method used to compute the non-divergent stream function 
is the procedure suggested by Thompson (1961), which begins 
with a field of isotachs and isogons. The Helmholtz 
equation is manipulated to form the equation: 

9 (Seegoy= v9 (8) 
By relaxing the field of the vorticity for suitable boun- 
dary conditions, the non-divergent stream function on an 
PsODaric SUrIace Is ODUalned. 

The numerical method and boundary conditions for 
solving this equation are those given by Hawkins and 
Rosenthal (1965). The bcundary values of the non-divergent 
stream function are made a function of the divergence of 
the wind, instead of the Zero condition imposed on the 
solution by the balance equation. The values depend on 
the derivative of the velocity potential (3%) (ny normal to 
the boundary, increasing outward). The boundary values of 
3% needed are calculated from an initial ~‘Y~field obtained 
from the solution of 

v*x =V°V, X=O on the boundary (9) 
where V:'V is the divergence calculated from the wind 
analysis. 

ine non-dliversenvuserreanm function calculaved by this 
method is what Baumhefner (1966) calls the kinematic stream 
function, and is the method that he felt was most accurate 


i cie Uropics. 


aa 


The numerical model accepts data at five levels (1000, 
ECOmNOOOWM4@0, andez00@embes) and utilizes the X, ¥)» P co= 
ordinate system. The basic inputs are, the observed tem- 
peratures, and the isogons and isotachs of the actual 
Wind ab Cache e—eeicmeasic grid congisus of 33 x P52 
points alone =the sZonaieemecradional and vertical directionse 
Actual grid values are inscribed at 27 grid points in the 
Zona lL GimeeGrone oul an arvbliilcial cyclical continuity es 
provided by extending the grid over 6 extra points by fit- 
ting a polynomial through the boundary conditions. This 
region is a numerical waste-box in the diagnostic studies 
but will become essential in the future prediction phase. 

A vertical staggering of the grid was used as recom- 
mended by Krishnamurti (1966a), and the values of YY and 
2 appear at 1000, 800, 600, 400 and 200 mbs while 6,7, 
oe and % = appear at 900, 700, 500, and 300 mbs. Other 
boundary conditions are, 0 =% = (%2%)20 at the northern and 
southern wall. Boundaries are eliminated in the Zonal 
direction by the use of the 6 additional points. W is taken 
tO be Zero at the surface and at the top of our atmospnener 
100 mbs. 

The wind velocity is used to obtain the vorticity, jam 
from the vorticity, the stream function is derived. The 
static stability is calculated from the temperature. These 
terms, together with a frictional effect from the strean- 
lines, will yleld the first approximation of the vertreau 


i 


Mou Lomi. 
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The omega equation, used for the first approximation, does 
not contain any terms that are a function of omega. The 
first guess of the vertical motion is then used to calculate 
the non-linear terms containing omega. The latent heat term 
is added at this point, since 1t is non-linear in nature. 

An iterative process is performed until the difference be- 
tween successive approximations is considered small. 

A brief recapitulation of the general form and effect 
of the main forcing functions is warranted. A complete 
discussion of all terms is included in Krishnamurti, 

Nogues, and Baumhefner (1966) and Krishnamurti (1966b). 

Basically, if “tne fercin=e function, F. 1s greater 
than zero, then in that region we can expect rising motions. 
There are exceptions to this rule, but generally this is true. 
Near the centers of maximum rising or sinking motions, FP, 
and GS. have been found to be inversevy proportional to 
each other (Krishnamurti, Nogues, and Baumhefner, 1966). 

Baumhefner (1966) has shown that the four main contri- 
ming forcime functions In the troplies are terms 1, 2, 5, 
and 6 (page 19). The first two terms are the simplest to 
interpret. The effect of differential vorticity advection 
(term 1) is that if in the upper air there is a region of 
strong positive vorticity advection by the non-divergent 
component of the wind, then Po ls greater than zero and 
rising motions would be found. The inverse holds for 
Sinking motions. The Laplacian of thermal advection 


(term 2) has the effect, that at centers where warm air 
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advection is a maximum by the non-divergent component of 
the wind, Py ls greater than Zero and again rising motions 
will be found. The inverse again holds. 

Mae other two terms, friction and latent heat, are vf 
great significance in the tropics and are explained in 
great detail by both Krishnamurti (1966b) and Baumhefner 
(1966). Generally, Pe Camiction)) Wiliecenmtmibutemr sama 
motions in regions of cyclonic relative vorticity and sink- 
ingumoti#ons, in regtonsyof anticyclonic. relative vortic uaa 
This effect will only be felt significantly at the 900 mb 
output level. A drag coefficient (C1) of 2.0 x 107? was 
used, and held constant in the horizontal. 

The latent heat term (term 5) contributes in the tno-= 
pics where the following conditions are met: 

1. The atmosphere is conditionally unstable. 
2. Net moisture convergence in vertical columns is 
greater than zero. 
3. Lhe retiaciveshumnaidity at. Lhegpanrtwenteikes aa 
greater than or equal to 60 percent. 
This term only effects levels above 900 mbs, as that level 
is used as the top of the boundary layer. It is apparent 
that the latent heat term will have its greatest contribu- 


tion to the rising motions where cumulus clouds are presene. 
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III, RESEARCH PROCEDURES AND PROBLEMS 

Initially, an attempt was made to produce a compre- 
hensive, five day analysis covering the major portion of 
the Pacific from 40°N to 40°S. Plotted, mercator, 
1:40,000,000 charts were obtained from the Hawaii 
Institute of Geophysics, University of Hawaii for 1 
through 5 March 1965. Map time is 00Z. This period 
was chosen to allow a comparison of our results with 
those of Smagorinsky, whose report will be published at 
a later date. Since hand analysis had to be employed to 
prepare the input data, time considerations became a 
limiting factor and case studies on 1 and 3 March were 
agreed upon. The model employed in this study contains 
27 GYid points in the zonal direction and 5 points»in ghe 
meridional direction spaced 2 degrees apart. This grid 
Was placed in the area of most dense data in the Western 
Pacific where the analysis was felt to be most reliable 
and a well defined, if rather large, ITCZ was observed 
on the nephanalysis. The grid is boundéd by the 128°E 
and 180° meridians and the 14°N and 14°S parallels. 
The area includes approximately 5.1 million square miles 
and contains sixteen reporting stations (Fig.2 ). 
Bleven of these are Radiosonde stations and 5 are Pibal 
stations, which report upper winds only. 

Streamline, isotach, and temperature analyses were 
prepared for the 1000, 850, 700, 500, 300, and 200 mb 


levels, as well as a surface pressure analysis. 
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Isogon values were extrapolated to grid points using 
a parallel motion protractor and tabulated to 1 degree 
accuracy. IJIsotachs were analyzed at 5 knot intervals and 
extrapolated to tenths of knots. Temperatures were ana- 
lyzed to whole degrees Centigrade and tabulated to 
tenths of degrees at grid points. Surface pressures were 
analyzed at one millibar intervals and tabulated to tenths 
of emillibars,. 

A preliminary moisture analysis was prepared but was 
not satisfactory for our purposes and is not included in 
this study. The reason for its omission was based on two 
considerations: 

1) The data was too sparce, and unreliable 
values were felt to have been generated 


at grid points between reporting stations. 
2) The nephanalysis revealed only the gross 
features of the moisture field. A subjec- 


tive scheme was considered, whereby typical 
relative humidities would be assigned to 
cloudy areas, and the moisture constructed 
from the humidities, however, lacking ac- 
tual cloud photographs this method was also 
considered likely to result in inadequate 
gradients in both the horizontal and vertical. 
The omission is a deficiency of this study since the 
Pelease of latent heat Is”a major contributor to the ver- 
mreal velocities in the tropics. The iniluence of vtervaim 


ts also neglected. 
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Large vertical velocities can be obtained through orograph- 
ic triggering mechanisms which would otherwise be unex- 
plained dynamically. 

Tiros IX negatives were obtained from NWRF Ashville, 
for the period 1 through 5 March 1965, but no useable 
frames Wevewobserveameovering the area of interest during 
the days analyzed. Nephanalyses were constructed from 
the Tiros IX catalog, and original facsimile nephs, 
scaled to the computer output grid size. Channel two, 

8-12 microns, Radiation Data was obtained from:Goddard 
Space flight Center and was analyzed for March 23rd, 

A minimum amount of subjectivity was employed in each 
of the analyses. The nephs were used as an aid in sparce 
data areas to locate major centers or regions of conver- 
gence. Wiederanders! Analysis of Monthly Mean Resultant 
Winds over the Pacific was also employed as a consistancy 
check. No reports were eliminated and gradients were 
extrapolated linearly in all fields. A few aircraft re- 
ports at non-standard levels were available and used if 
they could be extrapolated to the standard level without 
Significantly altering the basic analysis. Vertical space 
or time cross-sections to eliminate or smooth erroneous 
data were not constructed due to time and data limitations. 

The observed height fields were not analyzed for in- 
clusion in this report in a smooth form. Rough sketches of 
these were prepared by the authors to compare with the geo- 


potential heights constructed by the computer, on the basis 
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of input temperatures. The geopotential surfaces did not 
agree well with either the observed heights or the winds. 
Moreover, a continuous degradation was noted with increas- 
ing height. Baumhefner attempted to correct this in his 
study by subjectively modifying the geopotential heights 
to fit the observed winds, although weak temperature 
gradients may have been altered significantly. He remark- 
ed that using the observed, unadjusted geopotential heights 
to calculate thickness may lead to a more reliable calcul- 
ation of temperatures. The authors, however, did not 
attempt this experiment, although it appears worthy of 
investigation. 

To insure a smooth vertical temperature structure, 
successive averaging was accomplished, which resulted in 
the output temperatures varying slightly from input values. 

Since the machine accepts data at the even levels men- 
tioned previously, several techniques were employed to 
convert the data at the standard levels to the input 
levels within the computer. All data at the surface and 
200 mb levels were used without modification. The ten- 
perature at 850 mbs was reduced to the 500 mb value, using 
a tropical, standard atmosphere lapse rate. The 600 and 
400 mb levels were obtained using the mean of the 700 - 

500 and 500-300 mb levels, respectively. The same scheme 
was used for wind data except that the wind vectors were 
Separated into u and v components prior to averaging. The 


850 mb winds were used directly as the 800 mb values. 
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A vertical staggering of the grid is convenient due to 
the nature of the hydrostatic and W equations. The first 
part of the program, which does not include vertical veloci- 
ties, yields outputs consisting of, temperature, potential 
temperature and static stability at the 900, 700, 500, and 
300 mb levels. Geopotential heights, the stream function, 
and absolute vorticity appear at the 1000, 800, 600, 400, 
and 200 mb levels. The @& portion of the program provides 
vertical velocities at the 900, 700, 500, and 300 mb levels 
and we assume ™ = 0 at the surface and 100mbs which repre- 
sents the top of our atmosphere. The contributions to the 
total W by various forcing functions are obtained, which 
allow an examination of their effects in the ae The 
sensible heat term contributed nothing in this study since 
we assumed thermal equilibrium between the atmosphere and 
ocean in the tropics. A discussion of the output charts 
is made in Chapter IV. The program was solved on an IBM 
7094 computer and takes approximately one hour per case 
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IV. CASE STULIES 
0000Z 1 MARCH 1965 

synopt tes situation The most signifitant Teseures@in 
the surface pressure pattern on 1 March 1965, (Fig.3), are 
the low pressure center located in North Central Australia, 
and the trough extending from the Southern Hemisphere mid- 
latitudes to the equator. Two small high pressure centers 
are located in the Southern Hemisphere, one over New Guinea, 
and a second at 180° longitude. A secondary zonal trough is 
evident at 4°N extending from 145°E to 180°. The surface 
streamline pattern (Fig.5) agrees well with the pressure 
field. The winds are cyclonic around the Australian low 
With strong cyclonic shear and curvature, and the winds are 
anticyclonic around the high centered at 180°. The trough, 
which has been labeled Frontal on the nephanalysis, does 
not show up Glearlywin the wind field® bute this is,an area 
of very light winds and only a few reports Were available. 
An asymptote of convergence is located in the Zonal trough 
in the Northern Hemisphere. The winds to the north are 
extremely steady in direction but strong anticyclonic shear 
exists north of the isotach maximum at 6°N 164°EF. A cy- 
clonic eddy exists at 148°E just north of the equator. 
Equatorial Westerlies are formed as the N.E. Trades cross 
the equator and recurve west of 150°E under the influence 
Of a chamee inssien of tae Coriolis force. Ay swliwe ITCZ 
forms. The nephanalysis (Fig.4) depicts, in the eastern 


section (orbit 456), the clouds associated with the frontal 
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trough, an open area to the east and an ITC band located in 
the Northern Hemisphere in the area occupied by the zonal 
trough. The ITCZ is displaced far north of its mean posi- 
tion, which is at 8° South in February and early March. ihe 
western portion of the neph, orbits 457 and 458, shows most 
of the entire grid covered by clouds. The mostly open area, 
less than 5/10 coverage, in the northwest corner, is the 
result of a small high pressure ridge extending southward 
from higher latitudes. The ITCZ is an extremely wide band 
and is split into a northern and southern part. The south- 
ern band is not easily explained by the pressure field and 
we suspect a large terrain influence in this region. ‘The 
850 and 700 mb flow patterns (Fig. 6,7) are similar in most 
respects to the surface. At 500 mbs (Fig.9) the eddy which 
appeared north of the equator disappears and a large conver- 
gence area is found at 140°R 7°North. The Australian cyc- 
lone is still evident but the influence of the subtropical 
high pressure belts is becoming more apparent. At 300 and 
200 mbs (Figs.9,10) these belts are the dominating features. 
More variability is noted at 200 mbs than at 300 mbs. 

The observed temperature structure is very complex and 
appears to have little vertical continuity. Strong grad- 
ients may exist in localized areas which were not picked up 
in the analysis. The 850 and 200 mb temperatures are shown 
(Figs. 11,12) and they clearly show a reversal in the temper- 
ature gradient in the vicinity of the ITCZ. This feature was 


discovered by Schnapauff in the Atlantic and verified by 
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Fig. 5. Streamlines and isotachs at 1000 mbs, 00Z 1 March 
1965. Isotach interval, 5 knots. 





Fig. 6. Streamlines and isotachs at 850 mbs, 00Z 1 March 
1965. Isotach interval, 5 knots. 
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Fig. 7. Streamlines and isotachs at 700 mbs, 00Z 1 March 
1965. Isotach interval, 5 knots. 





Fig. 8. Streamlines and isotachs at 500 mbs, 00Z 1 March 
1965. Isotach interval, 5 knots. 
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Fig. 9. Streamlines and isotachs at 300 mbs, 00Z 1 March 
1965. Isotach interval, 5 knots. 





Fig. 10. Streamlines and isotachs at 200 mbs, 00Z 1 March 
1965. Isotach interval, 5 knots. 
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Fig. 11. Temperature analysis at 850 mbs, 00Z 1 March 1965. 
Temp. in degrees Centigrade. Isotherm interval, 
1 degree. 
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Fig. 12. Temperature analysis at 200 mbs, 00Z 1 March 1965. 
Temp. in degrees Centigrade. Isotherm interval, 
1 degree. 


38 


Simpson in the Pacific. Riehnl(1958)prepared a vertical 
cross-section of the typical temperature in the region of 
the equatorial trough, which is included for comparison 
(Fig.13a). Riehl also prepared a vertical cross-section 
of heights for this region (Fig.13b), and a comparison 
may be made with geopotential heights computed in the pro- 
gram from the surface pressure and vertical temperature 
structure. The 1000 mb geopotential height field (Fig.14) 
is shown as an example. 

An example of the absolute vorticity, computed from 
the observed wind field, is shown for the 1000 mb level 
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Fig. 14. Geopotential height at 1000 mbs, 00Z 1 March 1965. 
Heights in meters. Interval, 10 meters. 





Fig. 15. Aseilute vorticity at 1000 mbs, 00Z 1 March 1965. 
Interval 10 x 107° sec7!. 
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Results. The results obtained through the application 
Gaeunis model LO ou study are most encouraging, particu- 
larly in the Northern Hemisphere where better data cover- 
age Was available. The stream functions generated through 
the relaxation of the vorticity field at 800 and 200 mbs 
are shown for 1 March in Figures 16 and 17. A comparison 
was made of the non-divergent wind and the actual winds in 
a meridional cross-section. At 800 mbs, only the u con- 
ponent of tme actual wind is used since thevcress segewen 
passed through the circulation center at 148°E, and the 
directions were more variable at this level. At 200 mbs, 
the longitude was chosen so that the streamlines had 
thelipemesx imum = Zeomal values, and here, the Gotal isotacm 
value was used. These are plotted against the Zonal 
stream wind Uy and are shown in Figure 18. 

Baumhefner (1966) found that the non-divergent part 
of the wind was about 90% of the total wind and his results 
are clearly verified. Stream functions at the other levels 
yielded similar results. 

Dry adiabatic vertical velocities are shown at the 
900, 700, 500, and 300 mb levels, (Fig.19-22). An Wie 
was drawn in the areas of maximum rising motion on these 
charts such that it remained within the boundaries of the 
ITCZ indicated on the nephanalysis and at the same time 
Was vertically consistent. The Northern Hemispae mew. 


Was remarkable. 


42 





Jake 152 


Non-divergent stream functinn at 800 mbs, 
00Z 1 March 1965. Isoline interval, 100x107 


5 msec -. 
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Non-divergent stream function at 200 mbs, 
00Z 1 March 1965. Isoline interval, 10x1072 
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Fig. 18. A comparison of the input and non-divergent wind 
values. 
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Fig. 19. Total vertical velocity at 900 mbs, 00Z 1 March 
1965. Arrows indicate direction. Interval is 
10x1079 mbs secu, 
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Fig. 20. Same as above except 700 mbs. 
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Fig. 21. Total vertical velocity at 500 mbs, 00Z 1 March 
1965. -Arrows indicate direction, interval is 


10x107- mbs sec7!, 





Fig. 22. Same as above except 300 mbs. 


46 


Two points Were chosen in the vicinity of the BirCcZ, 
one at 164°E, 4°N in the eastern portion of the neph and 
a second at 138°E, 6°N in the western section, and a graph 
showing the variation of W in the vertical was constructed 
(Fig.23). Typical magnitudes were found on the order of 
fer icior teh ac lo 2 enie ceo frising) in the frictiong@i@layver, 
decreasing to less than 20 x 1072 mbs sec: above the 
frictional layer and with a secondary maximum in the mid 
and upper troposphere. This corresponds to Krishnamurti's 
theory that two maxima exist in the w field in the ver- 
tical, the lower caused primarily by the influence of 
frictional convergence, which drops off at intermediate 
levels, and the second maximum in the mid and higher trop- 


osphere, where more closed circulations are found. 
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Fig. 23. Vertical variation of W. 
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0000Z 3 March 1965 

Synoptic situation. By the third of March we find that 
the nephanalysis for 0251Z (Fig. 24), from Tiros IX orbits 
480-482, indicates that the Intertropical Convergence Zone 
is undergoing a major change with the northern section weak- 
ening and a second ITCZ intensifying to the south, so that 
there are actually two ITC's at 00Z on 3 March. The north- 
ernmost is defined mostly by the cloud picture and the use 
of the water vapor channel two data from Tiros VII indicates 
that it is most intense below about 20,000 feet. 

ime southern ITCZ rerwell defined by both elouds ame 
convergence of the wind field with the channel two data in- 
dicating it extends to greater than 30,000 feet. Since the 
Southern ITCZ is better defined by the wind? iweld, much of 
phe discussionof the stréamlines will cencern this area. 

It is interesting to note that the nephanalyses for 
the ist through the 4th of March show the continuous se- 
quence of events of the change in the ITCZ. This phenome- 
non has been observed before by Airborne Early Warning 


Squadron 4 which made daily flights in the Eastern Pacific 


equatorial region for a period of 13 months in 1957, during 
Which time it was observed that the ITCZ intermittently 
moves rapidly north, dissipates, and reforms in the south 
(Reth, 1958). 

The surface pressure analysis (Fig. 25) does not show 
the change in synoptic situation so dramatically, with the 


Omi Season front stvlll trailing down into the trepic region 
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Fig. 25. Surface pressure field, 00Z 3 March 1965. 
Interval, 1 millibar. 
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Fig. 26. Streamlines and isotachs at 1000 mbs, 00Z 3 March 
1965. streamlines, ----- isotachs. 
Isotach interval, 5 knots. 
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with low pressure extending northwest across the equator 
and a 1010 mb ridge located over the southern half of New 
Guinea. The streamlines and isotachs at the surface 

(Fig. 26) also support the pressure analysis, although we 
would expect a little more cross-isobar flow near the equa- 
tor because of the diminished Coriolis parameter. 

Tt is Ghe auvchors' opinion that the anticyclonic and 
cyclonic eddies to the south and north of New Guinea, re- 
spectively, are caused by the effect of the mountainous 
terrain of New Guinea on the wind flow. The cyclonic cir- 
culation in the north hemisphere has weakened considerably 
since the ist of March when it extended through several 
levels. On the 3rd it is only indicated by a slight wave 
at 850 mbs (Fig. 27) and disappears completely above that 
level. The southern ITCZ is well defined by the 850 mb 
streamlines and there is an excellent correlation with the- 
nephanalysis. 

The streamlines also indicate a reason for the weak- 
ening of the northern I17CZ,as divergence is indicate@ meas 
the wmicdiie of theszenme: 

As we go upward in the atmosphere (Fig. 28), we can 
see that all trace of the e€yclonic circulation’to the Miemel 
has disappeared but that the anticyclonic ridge over the 
seuth of New Guinea intensifies with height (Figs. §2ommueor 
and becomes a closed circulation by 300 mbs with outflow 
taking place. 


We also find that we now have some divergence taking 


place over the southern ITCZ. The Subtropical Highs have 
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Fig. 27. Streamlines and isotachs at 850 mbs, 00Z 3 March 
1965. streamlines, ----- tsotachs. 
Isotach interval, 5 knots. 
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Fig. 28. Streamlines and isotachs at 700 mbs, 00Z 3 March 
1965. streamlines, ----- tsotachs. 
Isotach interval, 5 knots. 
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Fig. 29. Streamlines and isotachs at 500 mbs, 00Z 3 March 
196 streamlines, ----- isotachs. 
Isotach interval, 5 knots. 








Fig. 30. Streamlines and isotachs at 300 mbs, 00Z 3 March 
1 streamlines, ----- isotachs. 
Rooracn interval, 5 knots. 





54 


moved equatorward and they can be seen clearly at 200 mbs 

(Fig. 31). The anticyclonic circulation persists at this 

level along with the divergence over the northern coast of 
New Guinea. There is now some convergence a little to the 
recibir ol, iia rea. 

By looking at the temperature structure at low, inter- 
mediate and high levels (Figs. 32-34), we can see the anti- 
cyclonic ridge appears cold near the surface and the cyclonic 
Gtreculation to the nortnhyis warm core. There is a pers#stent 
warm area to the southeast. What might be considered an 
anomaly is the fact that the region within 3 or 4 degrees 
latitude of the equator on both sides appears to be rela- 
tively cool at all levels compared to the surrounding area 
and this had a great effect on the calculated geopotential 
heights. When one notes that this is an area of consider- 
able cloudiness and moisture with high evaporation, it is 
more explainable (Riehl, 1954). The temperature at the 
levels shown are the smootned temperatures that were com- 
puted in the program and used in subsequent calculations. 

Output. Geopotential heights were built up from the 
surface pressure and the vertical temperature structure by 
bhne use of The standard avmosphere for the 3rd of March 
(Figs. 35-39). “W1lthough the heights are not used.in fur- 
ther calculations, they were done as a matter of interest 
to compare with Riehl's model (Fig. 13b) (Riehl and Malkus, 
1958). The geopotential heights do not indicate the ex-~ 


treme height anomaly that Riehl found at high levels over 
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Fig. 31. Streamlines and isotachs at 200 mbs, 00Z 3 March 
1965. Streamlines, —----— isotachs. 
Isotach interval, 5 knots. 
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Fig. 32. Temperature at 900 mbs, 00Z 3 March 1965. Temp. 
in degrees Kelvin. Isotherm interval, 1 degree. 
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Fig. 33. Temperature at 500 mbs, 00Z 3 March 1965. Temp. 
in degrees Kelvin. Isotherm interval, 1 degree. 
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Fig. 34. Temperature at 300 mbs, 00Z 3 March 1965. Temp. 
in degrees Kelvin. Isotherm interval, 1 degree. 
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Fig. 35. Geopotential height at 1000 mbs, 00Z 3 March 1965. 
Heights in meters. Interval, 10 meters. 
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Fig. 36. Geopotential height at 800 mbs, 00Z 3 March 1965. 
Heights in meters. Interval, 10 meters. 
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Fig. 37. Geopotential height at 600 mbs, 00Z 3 March 1965. 
Heights in meters. Interval, 10 meters. 
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Fig. 39. Geopotential height at 200 mbs, 00Z 3 March 1965. 
Heights in meters. Interval, 10 meters. 
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the ITCZ, instead low pressure was found. This is believed 
to be caused by the difference in the vertical temperature 
of Riehl's model (Fig. 13a) as compared to observed tempera- 
tures on 3 March. The warm anomalies up to 200 mbs were 

MOU found. 

A comparison of the 200 mb streamlines (Fig. 31) with 
the geopotential height at 200 mbs (Fig. 39) shows little 
resemblance to one another which might be an indication 
that an observed height analysis, if enough data were avail- 
able, would be considerably different from the geopotential 
Heomeno. sTnis is a further point to be taien into accoune 
when deciding whether to compute the non-divergent stream 
function from the geopotential height or alternately from 
tiemovserved Winds, as has been done here. 

etavic stability.  DUhesdry static stability outpu teat 
all levels indicatestwo main features, which are illustrated 
by the values at 900 and 300 mbs (Figs. 40, 41). One is 
that values are slightly positive everywhere within the 
wegen on une 3rd of Maren. The second feature indicates 
whaw we might expect in the tropics, that the gradient of 
dry static stability is very small. The same situation had 
been found on i March. “lhe cloud situation indicates that 
mao molsu Static stabiligyewould be a signifveant factor 
and that much of the area would be conditionally unstable, 
leading to large contributions to the upward vertical ve- 
iocrtCics trem the lavenuy meat forcing Tuncctcron. “hits would 


be mostly true in the areas of extensive cumulus cloud 


coverage, such as the ITCZ. 
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Static stability _at 1000 mbs, 00Z 3 March 1965. 
Interval, 10x10-3 meters2 mbs~2 sec-2, 
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Static stability at 300 mbs, 00Z 3 March 1965. 
Interval, 10x1073 meters2 mbs72 sec“, 
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Vorticity.s An example of the absolute Voc tici, sou. 
put using the actual winds is shown at 1000 mbs (Fig. 42). 
It must be remembered that the signs of cyclonic and anti- 
cyclonic vorticity are reversed in the Southern Hemisphere 
so that negative values indicate cyclonic vorticity south 
of the equator. This introduces some confusion, especially 
in areas such as the closed - 10 x lose Center au 178°R 
which is situated on the equator. The fact that the center 
extends into both hemispheres indicates anticyclonic vor- 
Hicity mm the northernweeet of the centerwand cyclonic im 
the southern half. This is not a contradiction if we take 
a look at the 1000 mb streamline chart (Fig. 26). The 
streamlines indicate a clockwise curvature and has a maxi- 
mum in the isotachs south of the easterly winds in that 
area. Both of these facts lead to just the picture that is 
mliiuserated by the vorttemyy center. Anticyclonic voremercy 
wemune north of the equavomwand cyecllonve vorticity to the 
SOU ul. 

The program takes the change of sign into account when 
Seompiieing vertical veloelules from vorticity effects sovemat 
the sign of © will be correct in both hemispheres. 

Stream function. Computed non-divergent stream func- 
tions are shown at 1000, 800, 600, 400 and 200 mbs (Figs. 
43-47) and a comparison with the streamlines at similar 
levels (Figs. 26-31) shows a striking similarity. A north- 
south cross-section was selected at random to see the actual 


Corre latwommen the Sra vacmcompared gro thatemound@on the Ist. 
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Fig. 42. Absolute vorticity at ee mbs, 00Z 3 March 1965. 
Interval, 10x 107” sec~ 
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Fig. 43. Non-divergent stream function at 1000 mbs, 
00Z 3 March 1965. Interval, 10x10-5 meters sec7l, 
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Fig. 44. Non-divergent stream function at 800 mbs, 
00Z 3 March 1965. Interval, 10x1072 meters sec™!, 





Fig. 45. Non-divergent stream function at 600 mbs, 
0OOZ 3 March 1965. Interval, 10x10-5 meters* sec7l, 
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Fig. 46. Non-divergent stream function at 400 mbs, 
00Z 3 March 1965. Interval, 10x10-5 meters® sec7! 





Fig. 47. Non-divergent stream function at 200 mbs, 1 
O0OZ 3 March 1965. Interval, 10x107? meters’ sec’. 
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The cross-section was taken at 140°E at 200 mbs. The actual 
winds were almost strictly Zonal in this area so the entire 
wind was taken and compared to the zonal non-divergent com- 
ponent calculated from Uy = - oa . The results (Ceiceee) 
show a remarkably close correlation, with some apparent op- 
posing winds from the divergent component between 3°N and 
ons and between 5° and 928%: Note that the true winds pass 
through a calm area about 8°S (Fig. 31) and change direction 
by almost 180° and the non-divergent component shows the 
same results. Near the centers of highs and lows we would 
expecye a Wargem contribution from the divergent partvorerne 
wind. 

Weruical velocity) malcmrEotal Verli calle veloc: ty cia 
Units of Ae S/o is shown at the output levels of 900, 
700, 500 and 300 mbs (Figs. 49-52) with the ITCZ's super- 
imposed as hatched double lines to observe the correlation, 
me any. 

It must be remembered that the vertical velocities 
shown here do not include sensible heat, which should give 
4 minor centribution in the tropics, and latent heat, which 
Has been sewn lo be a ma jo conikributor Gemihne upward ver— 
Pecalb Velocities 1n the trepics. We find a fairly good) ¢cer— 
relation with the cloud cover and the ITCZ's at 900 mbs 
(Fig. 49), with vertical velocities greater than 1 cm/sec 
upward and over 1.5 cm/sec downward. Rising vertical mo- 
tions are found over 95 percent of the analyzed ITCZ's with 


two intense centers in the north and a weaker one to the 
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Fig. 48. North-South cross-section of actual wind speed 
versus zonal component of the non-divergent wind at 140° E, 
200 mbs, 00z 3 Mar. 1965. 
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Fig. 49. Total vertical velocity at 900 mbs, 00Z 3 March 


1965. m—— sinking, -°-:+ neutral, ----- rising. 
ITCZ is cross-hatched. Interval, 20x1075 mbs sec”-. 
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Fig. 50. Total vertical velocity at 700 mbs, 00Z 3 March 
1965. sinking, --+--+ neutral, ----- rising. 
ITCZ is cross-hatched. Interval, 10x107> mbs sec7!, 
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Fig. 51. Total vertical velocity at 500 mbs, 00Z 3 March 





1965. sinking, -+++- neutral, -zz-- rising, 
ITCZ is cross-hatched. Interval, 10x10 mbs sec7l, 
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Fig. 52. Total vertical velocity at 300 mbs, 00Z 3 March 
1965. sinking, -°+*++ neutral, ----- rising, 
ITCZ is cross-hatched. Interval, 10x107 mbs sec -. 
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south, indicating that the northern ITCZ may be more active 
near the surface. The higher levels show a much poorer 
correlation with the ITCZ's. The 700 mb output (Fig. 50) 
does show rising motions over most of the southern ITCZ and 
part of the northern one. Maximum vertical velocities are 
of the order of .25 cm/sec. The correlation at 500 mbs 
(Fig. 51) is even worse except in one portion of the north- 
ern ITCZ where motions are as strong as .35 cm/sec, rising. 
The correlation at 300 mbs (Fig. 52) seems to be the worst 
of all, but the channel two data for the area indicates that 
clouds only extend above 30,000 feet in parts of the southern 
MfGZ. Vertical velocities are of the order of .6 cm/sec. 
The magnitude of the vertical velocities for 3 March at all 
levels are considerably less than for the ist. 

Latent heat. The authors feel that the poor correla- 
tion with the ITCZ's at levels above 900 mbs is due to the 
lack of the latent heat contribution. Baumhefner (1966) 
has shown that the contribution from the latent heat forcing 
function easily runs as high as .5 cm/sec and could contri- 
buve™mores These contributions to the upward motion would 
be found in areas with moisture convergence, R.H. > 60% and 
conditionally unstable conditions. These conditions would 
most likely be found in the areas of maximum cumulus cloud 
occurence which are right along our ITCZ's. Therefore, we 
feel that the addition of the latent heat team will bring 
the vertical velocities into a much closer correlation with 


the ITCZ's at levels above 900 mbs. 
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Inclusion of moisture for latent heat will be conducted 
this summer sin a follow wp study by Dr. T. N. Krishnamurti. 

Friction. Ihe higher correlation of the 900 mb level 
1s believed to be due in large part to the effects of fric- 
Plone! 12> ssaowemunemeontvribution fromethe friction 
forcing function at the 900 mb level and it can be seen that 
there is an excellent correlation with the ITCZ's and that 
friction is a large contributor to the total vertical mo- 
tion, running almost .4 cm/sec in some areas, which is over 
4O percent of the total velocity. The authors would agree 
with Baumhefnem that friction must be considered yon agpan 
with the forcing functions of the differential vorticity ad- 
vection by the non-divergent wind and the Laplacian of ther- 
mal advection by the non-divergent part of the wind as the 
NaniOngc oni abou toms. 

It might be mentioned that, initially, problems were 
encountered in computing friction at the 900 mb level and 
final friction contribution was not calculated untmibmeas 
the end of the research time period. As an alternative, in 
order to find the areas in which friction would have the 
most effect, a method to estimate the frictional vertical 
velocity at the topmetythe friction layer was utilized 
This method is discussed by Panofsky (1956) and concerns 
tne relationship of ~the frictional vertical veloolpge 


through the Ekman Theory, to Equation (10). 
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164 168 72 
Fig. 53. Frictional component of vertical velocity at 900 


mbs, 00Z 3 March 1965. m— sinking, ---+--+ neutral, 


— rising. ITCZ is cross-hatched. Interval, 
10x107? mbs sec7}. 





Fig. 54. Relative vorticity (for computation of friction) 


at 1000 mbs, 00Z 3 March 1965. —m—— anticyclonic 
OMe « ----- cyclonic vorticity. Interval, 
Sx107 sec7l, 
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locals acbie Vetmieal velocity producedmby friction 
aco the topwot the friction layer, f, is the surface stress, 
and $ denotes the average vorticity between the surface and 
poicwbCD Oh ticm: rlCLEOMaeayer. 

Panofsky shows that this can be expressed in the simpler 


fons 


. Mz @ 
Wi = a? Soy SIN 2X. 


Here, Ma is the eddy viscosity, eqs is the surface 


(1) 


geostrophic vorticity, and o&, 1s the angle between the sur- 
face wind and surface isobars. The quantity /42 can also be 
eliminated by means of the relation between the height of 
the friction layer, Af , and Ma . The result is equation 
CUZ 


S ogs Sin AX. , (128) 





We = YT+2%o 

An estimate can be calculated for the frictional ver- 
tical velocity by subtracting the Coriolis parameter (f) 
from the absolute vorticity to obtain the relative vorvleiey 
(Fig. 54). This then can be used in Equation (12), along 
with a constant angle, &,, over the ocean. 

An example would be 

for $230 X10 /Ssec, HE IKM , AND xo 350", 
wy 2 0:35 am /sec 

A comparison of Fig. 53 with Fig. 54 shows a high de- 
eres of correlation, if One remembers that negative value. 
of 3 in the Southern Hemisphere represent cyclonic vortici- 


ty and we would have rising motions where they are located. 
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Other terms. Baumhefner (1966) found that in many in- 
stances the contributions of the vorticity advecviom ang 
the Laplacian of thermal advection were of opposite sign. 
The same results were noted here, as illustrated by a conm- 
parison of the individual partial vertical velocities due 
to these forcing functions at 700 mbs (Figs. 55, 56). 

The fact that the first two forcing functions, men- 
tioned above, are generally the only major contributors to 
the vertical velocities above 900 mbs, in the absence of 
latent heat, is illustrated by a closer look at the 700 mb 
level. Fig. 57 shows the results if the sum of the contri- 
butions from vorticity advection and the Laplacian of thermal 
advection is subtracted from the total vertical velocity. 
It can be seen that the contribution of the sum of the re- 
maining eight forcing functions that were calculated is al- 
most negligible. 

Mertical stmucture. 4 point along the northern ITCzZ 
was selected for the 3rd of March, as was done for the ist, 
to check the vertical structure of the magnitudes of the 
vertical velocities. The results at 8°N 144°E are shown in 
Fig. 58 and show the same basic structure as on the ist. 
There are large values near the surface, due to friction, 
overlaid by small values at 700 mbs and larger values at 
meanem Levels before dropping off to zero at the top ef wae 


atmosphere. 
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Fig. 55. Vertical velocity component due to differential 
vorticity advection by the non-divergent part of 
the wind, 700 mbs, 00Z 3 March 1965. sania 
eeee* neutral. Interval, 10x1075 mbs sec™ 
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160 164 168 


Vertical velocity component due to the Laplacian 
of thermal advection by the non-divergent part of 
the wind, 700 mbs, 00Z 3 March 1965. sinking, 
"ces; neutral, ----- rising. Interval, 10x10-5 mbs 





sec 





Net contribution of vertical velocity from eight 
forcing functions after subtracting vorticity ad- 
vection and Laplacian of thermal advection compo- 
nents from the total vertical velocity at 700 mbs, 
OOZ 3 March 1965. Sinking, °**** neutral. 
Interval, 10x1075 mbs sec-!, 
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Fig. 58. Pressure versus vertical velocity (w) at 8°N 
144°E, 002 3 March 1965. 
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V. CONCLUSIONS AND RECOMMENDATIONS FOR FURTHER STUDY 

Several conclusions and recommendations are made based 
on the results of the study. 

The relation between the calculated geopotential height 
fields and the stream fields becomes increasingly degraded 
with height, with the resemblance becoming negligible by 
200 mbs. It appears that computation of the actual geopoten- 
tial height field would be far superior in representative- 
ness if this field were to be used in further calculations. 

The close correlation of the non-divergent part of the 
wind to the actual wind illustrates the validity of current 
barotropic progs currently used by the weather bureau in 
Hawaii and by ESSA. These progs use just the non-divergent 
part of the Mee oe ENG mobs, from 30°N to 30°S. The non- 
divergent winds are obtained from the actual isogons and 
isotachs with calculations leading to Y*W. The resulting 
Winds are just moved ahead by prognosis. The results ob- 
tained by this method, since the non-divergent part of the 
wind represents about 90 percent of the actual wind, is 
probably better than what can be done by other methods until 
more complete information on winds and cloud coverage is 
available. 

It is felt that the lack of a moisture input consider- 
ably degraded the correlations that might be expected be- 
tween the vertical velocities and the ITCZ at levels above 
900 mbs. Previous studies have shown that the latent heat 


term is a @eonsiderable contributer to the upward vertical 


(ee 


MOMiCot@——e Sieermally In resions of extensive cumulus activity, 
such as the ITCZ. It is felt that the rising motion centers 
would have been more centrally located near the ITCZ's had 
Picea tenbeieay scorn been calculated.. Further studies 
should definitely include this term in the tropics. 

Iwo aids to the few dew point temperatures available 
Peorycompuving Molsvumesare cloud pictures and the channel 
two water vapor data. Areas and types of cloud coverage may 
be obtained from cloud pictures and the channel two data is 
a good indication as to what level the clouds extend, thus 
being able to drop off the moisture content above that level. 
Krishnamurti has suggested an exponential decrease in mois- 
ture with height and this will be used in future studies as 
an alternate method of obtaining moisture. Nephanalyses 
appear to be too gross for use in numerical calculat Yonseadgs 
this 1s especially true in the tropics. 

An illustration of the difficulties with the nephanal- 
yses is the case of 1 March 1965. Fully 90 percent of the 
grid is covered by five-tenths or more clouds. Perhaps a 
more reliable interpretation of results can be obtained by 
selecting a day when photographs show a better defined, 
single ITCZ closer to its mean position for a particular 
season. 

Terrain effects, although not included in* this stu, 
should be added in any region where significant land eleva- 
tions occur. The effect of the high mountainous terrain of 
New Guinea could have made a significant difference in the 
final analysis. 
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The preliminary results of the vertical variation of 
@ at several points indicate that the tropical atmosphere 
may be characterized by large values of vertical motion 
near the surface, due to friction, small values at interme- 
diate levels with much non-divergent flow and large values 
again at some higher level where the effect of a greater 
number of closed circulations is realized. The vertical 
motions are assumed Zero at 100 mbs. This would indicate a 
vertical picture of strong inflow into a system near the 
surface with weak inflow or perhaps outflow near 700 mbs, 
inflow again above that level, and outflow near the top of 
the atmosphere. Further investigation into this subject is 
recommended. 

It should be mentioned that a short range forecast is 
being made for this area in the future, utilizing the re- 
sults obtained here in a primitive equation prediction 
model. This is being done in a manner correct for small 
Rossby number by Dr. T. N. Krishnamurti. 

Further observational work, both through the use of 
conventional data and through the use of satellite cloud 
photography and infrared radiometric probes will undoubtedly 
provide better grounds for parameterization of the sub-grid 
scale phenomenon, and continue to give us a better under- 


standing of the tropical atmosphere. 
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